A two-way algorithm for the entanglement problem 
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We propose an algorithm which proves a given bipartite quantum state to be separable in a finite 
number of steps. Our approach is based on the search for a decomposition via a countable subset of 
product states, which is dense within all product states. Performing our algorithm simultaneously 
with the algorithm by Doherty, Parrilo and Spedalieri (which proves a quantum state to be entangled 
in a finite number of steps) leads to a two-way algorithm that terminates for any input state. Only 
for a set of arbitrary small measure near the border between separable and entangled states the 
result is inconclusive. 
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The question of whether a given quantum state is en- 
tangled or separable is both of fundamental interest, and 
of relevance for the implementation of quantum informa- 
tion processing tasks. The separability problem has stim- 
ulated many ideas for partial solutions: A sufficient con- 
dition for separabilityis given by the vicinity of the state 
to the identity yj, |2, |j| • A necessary condition for sepa- 
rability of a given state is that it fulfills the criterion of 
the positive partial transpose (PPT) [4|. Entanglement 
witnesses provide sufficient criteria for entanglement [5J. 
However, the separability problem has been shown in |6( 
to be in the complexity class NP-hard, and no complete 
solution is known yet. An improved algorithm for the 
separability problem, based on entanglement witnesses, 
was recently proposed in 0- In this Letter we suggest 
an algorithm that extends and complements the recent 
algorithm by Doherty, Parrilo and Spedalieri 

The separability problem is defined as follows. A quan- 
tum state p which acts on a bipartite, finite-dimensional 
Hilbert space Ha <8> Hb is separable iff there exists a 
set of pure product states |e<)(e<| (8 |/»)(/»|, and a set 
of real positive numbers pi with YliPi — 1, such that 
P = J2iPi\ e i)( e i\ ® I /»>(/»! 9J- Tuis property can be 
reformulated such that p has to lie within the convex 
hull of some pure product states. Furthermore, it is 
known that a separable p is in the convex hull of at most 
L := (dim Ha dim Hb) 2 pure product states 0. So 
it remains to show whether for a given p there exist L 
(not necessarily pairwise different) pure product states, 
such that p is in their convex hull. However, searching 
for these L states in the set of all quantum states would 
mean to search through a set of infinitely many uncount- 
able states. 

One of the most advanced solutions to the separability 
problem was recently introduced by Doherty, Parrilo and 
Spedalieri in They presented an iterative algorithm 
(denoted as A\ in the remainder of this text), which is 
based on symmetric extensions of a given quantum state, 
such that this algorithm terminates after a finite number 
of iterations iff the state is entangled. However, if the 
state is separable, the algorithm A\ does not terminate. 



From an algorithm-theoretic point of view this is not sat- 
isfactory: not having terminated after a finite time does 
not yield any information about the properties of the 
state. 

We suggest an algorithm which provides a solution for 
this problem and closes the gap in the above algorithm, 
because it detects a given separable state after a finite 
number of steps. Applied in parallel with the algorithm 
in ||, the combined algorithm then terminates after a 
finite time: one of the two tests certainly terminates, as 
every state is either entangled or separable. Our main 
idea is that it is sufficient to restrict ourselves to a count- 
able subset of pure states, rather than searching through 
all (uncountable) pure states. 

Let us start by providing the mathematical back- 
ground. For each separable state p there exists by def- 
inition a decomposition p = J2i=i Pi\ e i)( e i\ ® !/>){/> 1 1 
with pi > and YliPi = 1- A central idea of our 
approach is that the probabilities pi are not needed 
for deciding whether a state is separable or not - only 
the pure states {|ei),|/i)} in the decomposition are es- 
sential. Thus we can rephrase the separability defini- 
tion in the following way: A state p is separable, iff 
there exists a set of projectors onto separable pure states 
c := {|ei)(ei|®|/i)(/i|, \e L )(e L \(g>\f L )(f L \}, such that 
p lies in the convex hull of the elements of c. However, a 
straightforward search through all sets containing L pure 
product states - while easily parametrized - is impossible, 
as there are not only infinitely many such sets, but they 
are even uncountable. 

We will show in the following that it is (in the generic 
case) enough to restrict oneself to a countable subset C of 
all L-tuples of pure product states, which is dense within 
all pure product states. Within a countable set C there 
exists by definition a sequence {ci} in C that covers C 
completely. Thus it is possible to formulate an iterative 
algorithm that passes all L-tuples of pure product states 
in C in the limit of infinitely many steps. Furthermore 
one can use the well-known feature that for every element 
c in C there exists a finite number i such that the value 
Ci at the i-th step of the sequence equals c. Therefore 
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for a given (generic) separable p the iterative algorithm: 
Check whether p is in the convex hull of the L-tuple Ci 
terminates after a finite time. 

Let us put the mentioned ideas on mathematical 
grounds and first prove the statement that it is suffi- 
cient to search within a dense subset. Although some of 
the arguments given below hold for general convex sets, 
we will restrict ourselves here and in the following to 
sets of operators ( "states" ) acting on a finite-dimensional 
Hilbert space, where we use the Hilbert-Schmidt norm. 
The distance between two vectors from the Hilbert space 
is given by d(\tf>), \<f>)) = y/M - WMEMh and be " 
tween two operators d(a, b) = yTr[(at — Wj(a — b)]. 

Def. 1: A subset B C A is called dense in A if every 
a 6 A can be written as the limes of a sequence {b n } in 
B. 

Let us summarize some facts about convex sets: 

1. A set A is convex if for any finite number I, any 
ax,...,ai € A, and all A; > with Ei^j = 1, the sum 

Ei=l ^i a i € X. 

2. The convex hull of a set X is the smallest convex set 
that contains X. The convex hull of the set X will be 
denoted as conv A. 

3. For a convex set X we will denote the border SX as 
all points a € A for which 3b G X such that for all rj > 
the point (1 + r/)a — r\b does not belong to X. 

Lemma 1: Let A = c\A (where c\A denotes the clo- 
sure of A) be the set of extremal points of a convex set A, 
and B C A be a dense subset of these extremal points, 
then the convex hull of B is dense within the convex hull 
of A. 

Proof: We want to prove that conv^4= clconvB, so 
we have to show both inclusions. 

"c": Suppose one has x S conv^4 then there exists by 
definition a set of elements a,i in A and Xi > with 
= 1' such that x = El=i^ a 2- Since B is dense 
in A there exists for every Oi a sequence {bij} in B such 
that the limit j i— * oo of this sequence is a^. Due to the 
additivity of the limes x is in the closure of conv_B. 
"D" : If A is closed then conv A is also closed, and con- 
tains conv B. But cl conv B is the smallest closed set that 
contains conv-B, so clconv-B C conv A. □ 

Lemma 2: Given a convex set X within the set of 
states, with the finite set of extremal points ai,a2, ...a/ 
and a point x that is in the interior of X. Then there is 
an e > such that for all points a' with d(ai, a') < e the 
point x is also in the interior of the convex hull of the 
points a', a%, ...a/. 

Proof: A "face" of a convex set with dimension D 
is a subset of the border, that is the convex hull of D 
affinely independent elements. Since x is in the interior 
of convjai, ..a;}, there is an e > such that all points 
y in A for which d(x, y) < e are not on the border of 
A. In other words: x is not in an e-surrounding of the 
border. We note now that changing a\ to a' only affects 




FIG. 1: Constructing a "generalized diamond", i.e. a 
regular crosspolytope. 

those faces that have a\ as one of the affinely independent 
extremal points. When a' is taken such that d{a\, a 1 ) < e, 
the involved faces will remain in an e-surrounding of the 
faces of the original set of extremal points. Therefore the 
point x remains in the interior of the new convex set. □ 
Theorem 1: Given two convex sets B C A, such that 
B is dense within A, then all points in the interior of A 
are in the interior of B. 

Proof: Consider an arbitrary point a in the interior of 
A. Then Be > s.t. all points with a distance smaller 
than 2e from a are in A. Let us choose 21 points defined 
by Xi = a + e/3j, xi + i = a — efli, with i = 1, I, where 
{/3i, f3i} is an orthonormal basis in our space. Each of 
these points has the distance e from a, and is therefore 
in A. The convex hull of these 21 extremal points is a 
"generalized diamond" , i.e. a regular crosspolytope (see 
Fig. (IJ, and a has the distance y/le > from any face of 
this diamond. Using Lemma 2, we can shift one extremal 
point of the diamond after the other, such that the shifted 
points are in B (as B is dense in A), while keeping a in 
the interior of the new diamond. As all new extremal 
points belong to B, a is in the interior of B. □ 

Let us apply these general properties of convex sets 
to the problem of proving the separability of a given p. 
We will parametrize a countable subset of pure product 
vectors that is dense within all pure product vectors, and 
show that this leads to a dense subset of corresponding 1- 
dimensional projectors. Here, we have to distinguish the 
cases where p has full rank, or does not have full rank. 

We first study the case that p has full rank. The 
set of all pure product vectors is parameterized by fix- 
ing an orthonormal basis {|1)a, — , \ u )a} in Ha, where 
dim Ha = n, and an orthonormal basis {|l)s, 
in H.b, where dim7is = m. Then every pure product 
state can be written as 1^) = |a) ® \b) = (YliXi\i)A) ® 

(E,-Mili>fl)i wher e E* l^| 2 = 1 and Ej N 2 = 1- So 
the pure separable states are parametrized by the set 

G= {Ai,...,A„,Mi,...,Mm|EiM 2 = 1,E, M 2 = 1} of 
n + m complex coefficients with the two normalization 
constraints. 
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We now restrict the coefficients to complex quantities 
expressed with rational numbers Q. Note that the ratio- 
nal numbers are a countable set and dense in the set of 
real numbers R. However, due to the normalization con- 
straint we cannot simply consider the subset of G where 
all coefficients are of the form | + if (where p, q, r and 
s are natural numbers). We solve this problem by em- 
bedding the normalization constraints explicitly, choos- 
ing the subset 
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The subset G is dense within G, since for every element 
g = (A' l7 ...,Mm) m G there is an element g = (Ai, ...,p m ) 
in G that is arbitrary close to <?, when the distance 
is defined as d(g,g) = y/(X[ - A x ) 2 + ... + (fi' m - p m ) 2 . 
Furthermore G is countable, since it is a subset of 

QX2(n+m-l) 

Obviously the product vectors parametrized by G are 
dense within all product vectors in the Hilbert space, 
since the distance d(g, g) is equal to the distance induced 
by the Hilbert-Schmidt norm. 

Lemma 3: If a sequence of normalized vectors \ipi) 
converges towards \<fi), then the corresponding projectors 
\4>i)(4'i\ converge towards the projector \4>){(j)\- 

Proof: The distance betw een \tpi) and \<fr) is 
d(\i>i),\(f>}) = y/2(l - 3?[(^|^)]), where we denote by 
3?[('i/'i|(/>)] the real part of the scalar product. Since the 
\ipi) are converging towards \4>), for every e > there 
exists an i$ such that for all i > iq the distance d < e. 
This implies that 9fJ[(V> l |</>)] > 1 - e 2 /2. The distance of 
the corresponding operators is calculated as 

dmiuwm = v^-m\m (i) 

< v/2(l-(3?[(^|^)]) 2 ) < e^/2~7/2 . 

Thus \il>i)(il>i\ converges towards \<fi)((t>\. □ 
We now study the case that p does not have full rank. 
Since the states with lower rank form the border of all 
states, they will not necessarily be in the convex hull of 
the previously defined countable set. Thus we have to 
define the set G in a different way. 

Let r be the rank of p. We restrict ourselves to the 
(r 2 — l)-dimensional space of Hermitean operators with 
trace 1 that are supported at most on the range of p. 
Therefore the maximal number of extremal points needed 
to find a separable decomposition is given by L := r 2 . We 
know that p is in the interior of the space spanned by the 
projectors whose corresponding vectors are in the range 



of p. Thus, in the case of less than maximal rank we do 
not have to check whether p is in the convex hull of all 
separable pure states, but whether p is in the convex hull 
of all product projectors, whose vectors are in its range. 
Note that the existence of "enough" product vectors in 
the range is a necessary, but not sufficient criterion for 
separability, since any state of full rank has all pure states 
in its range, independently of its separability property. 

Given a state p and its spectral decomposition p = 
J2i=i Pi\4>i) {4>i\, a vector \ip) in the range of p can be 
written as = Ym=\ Ai|<fc)j with complex coefficients 
Ai and Y^i=i |Ai| 2 = L A pure bi-partite state 1^) is 
separable iff Tr J 4(TrB|V , )(V'l) 2 = 1- Therefore the coef- 
ficients for all pure product states in the range of p are 
the roots of a polynom of fourth order. These roots can 
be obtained numerically. The conditions for the A^ are 
summarised as follows: 
1. A, e C o A, t = e 27rie '|Ai| 

2- E[=i |A,| 2 = 1 |A,| = ^-EUilA,! 2 
3. Tr A (Tr s |V)(V|) 2 -l = 0, i.e. 

n m r r 

E E E E hK^i'K'i3AiB\<t>i){<j>k\3' A iB) ■ 

j,j'=l i,*'=l (,fe=l l',k'=l 

■OVsl^X^liA4)-i = o (2) 

We parametrize the product vectors in the range of p 
by 2r — 2 real parameters. Now we once again construct 
a dense subset of these parameters by choosing rationals 
in the form Aj = ^ L e 27rl ~ , where pj, q^, Ti, Sj G N and 
Pi < qi.Ti < Si. As previously, this subset is dense in all 
product vectors and therefore dense in the correspond- 
ing projectors. Furthermore it is countable, due to the 
countability of the rationals. 

Having found a countable subset that is dense in all 
product states, this immediately leads to a countable sub- 
set G of L-tuples of product states: each fixed L-tuple q 
can be provided with a "finite address" i. We can restrict 
ourselves to tuples with affincly independent elements: 
An affincly dependent set can be reduced to an affincly 
independent subset. For a decomposition with less than 
L product states, one can extend the corresponding tuple 
to an L-tuple by adding affincly independent entries. Ob- 
viously the state is still in the convex hull of the extended 
tuple, and will be detected as separable at its "address" . 

Our arguments lead to an algorithm for the detection 
of a separable state, in the following denoted as A^- One 
walks step by step through the countable set of L-tuples 
G. The ith-clement of G is c 4 = {r[ i] , r| l) }. One 

checks if the {t[ , t^} are affinely independent, and if 
they are not, one moves to the next element, i.e. i 
If the elements of are independent, one checks whether 
p belongs to conv{r}*\ r|^}. 

The check whether p is in conv{r}*\ r^} is per- 
formed as follows: one chooses L — 1 different elements 
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out of the L given ones, and finds the normal to the 
hyperplane defined by these elements. The state p is on 
the "same side" of the hyperplane as the remaining point, 
if sign(Tr[£ n p\) = sign(Tr[£„(T r - T h )]), where r r is the 
remaining point, and t/j is a point in the hyperplane. If 
the two signs are different (i.e. p and the remaining el- 
ement are "on different sides"), then p does not belong 
to coiw{t[ 1 \ ...,t^}. This test is performed L times for 
all possible choices of L — 1 elements from a. If p is in 
each case on the same side as the remaining point, then 
p is in conv{r} 4 \ ...t^ }, and therefore separable. In this 
case the algorithm terminates. Otherwise one continues 
with the next step, i.e. i i— ► i + 1. 

Combining the two algorithms A\ and A2 by running 
them parallel (in an iterative way) is already a big im- 
provement over Ai, since the combined algorithm ter- 
minates after a finite time for all input states that are 
entangled and all input states that are in the interior 
of the separable states. This leaves a set at the border 
between separable and entangled states, where the com- 
bined algorithm cannot be trusted to terminate after a 
finite time. This problem can be solved as follows. 

A state p on the border between separable and en- 
tangled states has the property that for all < 77 < 1 
the operator p e = (1 + rf)p — 77 1 does not belong to 
the separable states Due to convexity the opera- 

tor p s = (1 — rj)p + 77 1 is separable. If for all 77 > the 
operator p e is non-positive, then p is not of full rank - 
a case that we already studied above. Thus, the only 
possibility for a state to be on the border between sepa- 
rable and entangled states is that 3 770 > s.t. V77 < 770 
the state p e is positive. Note that until today there is 
no algorithm known for the decision whether a state is 
on the border between separable and entangled states: If 
this border would be known completely the separability 
problem would be solved. 

The above property can be used for closing the 
termination-gap in the combined algorithm described 
above, by extending it in the following way to the fi- 
nal algorithm A: Take some small, but fixed 77 > 0, such 
that (1 + ij)p — 77 1 is a positive operator. Then set two 
flags fi, fa to FALSE. These are global flags and are not 
changed at any step of the algorithm, unless mentioned 
explicitly. In the i-th step of the algorithm: 

1. do the i-th step of A\ for p, 

2. do the i-th step of A2 for p, 

3. do the i-th step of Ai for the state (1 + r\)p — 77 1, 

4. do the i-th step of Ai for the state (1 — rj)p + 77 1. 

If Ai detects {l+rj)p-r) 1 in 3., set fa to TRUE (from this 
point on it will stay TRUE). If A2 detects (1 — r\)p + 77 1 
in 4., set fa to TRUE (from this point on it will stay 
TRUE). 

The termination criteria for A are given as: 
a) If the A\ test detects p in 1., then p is entangled and 
A terminates. 



b) If the A2 test detects p in 2., then p is separable and 
A terminates. 

c) If both fa and fa are TRUE, then the state is in the 
77-surrounding of the border between separable and en- 
tangled states, and A terminates with this information. 
Otherwise do the step i 1— ► i + 1. 

This algorithm terminates after a finite time for any 
initial state. Here, the outcome c) does not give any 
information about the state being separable or entangled, 
but just the knowledge that the state is "close" to the 
border. However, we point out that the surrounding of 
the border which leads to an inconclusive outcome can, 
in principle, be made arbitrarily small. 

In summary, we have presented an algorithm for the 
separability problem, which complements the algorithm 
of Doherty, Parrilo, and Spedalieri. Their algorithm de- 
tects an entangled state after a finite number of steps, 
but does not terminate for separable states. Our algo- 
rithm, on the other hand, detects a separable state after 
a finite number of steps, but does not terminate for an 
entangled state. The connection of the two algorithms 
terminates for all input states. In the case of the ini- 
tial state being close to the border between separable 
and entangled states, our algorithm terminates with an 
inconclusive output. 
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